NASA Contractor Report 3069 



loan copy-'^MM 
K! »tla No 


Feasibility of Combining 
Linear Theory and Impact Theory 
Methods for the Analysis and 
Design of High Speed Configurations 



D. Brooke and D. V. Vondrasek 


CONTRACT NAS 1-1 5074 
DECEMBER 1978 




TECH LIBRARY KAFB, NM 



NASA Contractor Report 3069 


TECH LIBRARY KAFB, NM 

11111111 

□□blBflb 


Feasibility of Combining 
Linear Theory and Impact Theory 
Methods for the Analysis and 
Design of High Speed Configurations 


D. Brooke and D. V. Vondrasek 
McDonnell Douglas Corporation 
St. Louis, Missouri 


Prepared for 

Langley Research Center 

under Contract NAS1-15074 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Office 


1979 




TABLE OF CONTENTS 

Section Page 

SUMMARY 1 

INTRODUCTION 1 

ANALYSIS 2 

Shortcomings of Linear Theory and Impact Theory in 

Mach 4 to 8 Range 2 

Linear Theory 2 

Impact Theory 3 

Combined Impact Theory/Linear Theory Approach ... 4 

VALIDATION OF COMBINED IMPACT THEORY/LINEAR THEORY METHOD 5 

70° Swept Flat Plate at Mach 6.0 5 

76° Swept Wing Comparisons 7 

Additional Comparisons on the 76° Swept Wing .... 13 

Ames All-Body Models 20 

Closed Form Solution for Lift, Drag and Moment ... 22 

INVERSE DESIGN PROCEDURE 25 

CONCLUDING REMARKS 28 

APPENDIX 29 

REFERENCES 3 7 


iii 



Figure No. 
1 


LIST OF ILLUSTRATIONS 


Page 


Comparison of Region of Influence as Defined 
by Mach Angle (Linear Theory) and Shock Angle 
on 10° Cone 3 

2 Pressure Comparisons on 70° Swept Flat Plate, 

Mach 6.0,a=8° 6 

3 Effect of Scale Factor on Gentry /Woodward 
Solution for 70° Swept Flat Plate, Mach 6.0, 

cc = 8 ° 7 

4 Pressure Distributions on 76° Swept Wing (Mach 

4.6) 8 

5 Pressure Differences on 76° Swept VJing Using 

Current Method (Mach 4.6) 11 

6 Pressure Distribution on 76° Sweet Wing (Mach 

3.5) 13 

7 Pressure Distributions on 76° Swept Wing (Mach 

2.3) 16 

8 Comparison of Linear Theory, Impact Theory and 

Combined Approach for Pressures on 76° Swept Wing 
(Mach 2.3) 19 

9 Lift and Moment on 76° Swept VJing at Mach 2.3. . . 20 

10 NASA AMES All-Body Models Bp and B 5 21 

11 Lift and Moment Predictions on Ames All-Body 

Model B 5 , Mach 2.0 21 

12 Lift and Moment Predictions on Ames All-Body 

Model B 5 , Mach 5.37 22 

13 Closed Form Expressions For Lift, Moment and 

Drag-Due-to-Lift For Uncambered Configurations . . 23 

14 Closed Form Predictions - Ames Bodies Bp and B 5 , 

Mach 5.37 24 

15 Closed Form Predictions - Ames Bodies Bp and B 5 , 

Mach 7.38 24 

16 Optimization Matrix For Design Procedure 26 

17 Steps in Design Procedure 27 


IV 



SUMMARY 


This investigation examined the feasibility of combining 
elements of linear theory and impact theory to provide improved 
aerodynamic predictions for the analysis and design of high 
speed configurations in the Mach 4 to 8 range. Specifically, 
the aerodynamic influence coefficients calculated using an exist- 
ing linear theory program were used to modify the pressures cal- 
culated using impact theory. This combined approach was used to 
calculate local pressures and loadings over several wing-alone 
configurations. Comparisons with experimental pressure data show 
that the combined approach gives improved predictions over either 
linear theory alone or impact theory alone. The results of the 
study show that the combination of linear theory and impact 
theory not only removes most of the shortcomings of the indivi- 
dual methods as applied in the Mach 4 to 8 range, but also pro- 
vides the basis for an inverse design procedure applicable to 
high speed configurations. 


INTRODUCTION 

Impact theory and linear theory have provided the basic 
aerodynamic analysis tools used by the vehicle designer in the 
supersonic and hypersonic speed ranges. The linear theory 
methods, e.g.. References 1 to 3, apply in the supersonic range, 
while impact theory methods. Reference 4, are used in the hyper- 
sonic range. Strictly speaking, neither method is appropriate 
in the Mach 4 to 8 range, and both methods may exhibit anomalous, 
and often severe, failures when applied in that speed range. The 
following sections identify the specific nature of these fail- 
ures, and show how linear theory and impact theory can be com- 
bined to provide improved aerodynamic predictions in the Mach 4 
to 8 range. Comparisons with experimental pressure data are 
presented to show that the combined linear theory/impact theory 
approach gives significantly improved results over either theory 
applied by itself. 

In Reference 5, linear theory was applied in the inverse 
design mode to investigate the potential gains due to wing 
camber at high Mach numbers. The results indicated little or 
no payoff above Mach 5. The present analysis indicates that 
unmodified linear theory applied in the inverse design mode 
will underpredict the effectiveness of wing camber in reducing 
drag-due-to-lif t in the Mach 4 to 8 range. In this investiga- 
tion, the optimization equations are rederived to properly account 
for the coupling of the thickness and lifting terms, which is 
excluded from the linear theory analysis. In addition, a pro- 
cedure for the application of the combined linear theory/ impact 
theory approach in the inverse design mode is illustrated. 



In this study, the original Ames version of the Woodward 
program. Reference 1, and the Hypersonic Arbitrary Body Program 
or Gentry program. Reference 4, were modified and used for all 
basic calculations. The Woodward program provides the linear 
theory influence coefficients, which are used to 'correct' the 
two-dimensional, non-interference impact pressure distributions 
calculated using the Gentry program. The modified pressure dis- 
tributions are then integrated over the Gentry geometry model 
for aerodynamic forces and moments. 

ANALYSIS 


Shortcomings of Linear Theory and Impact Theory 
in Mach 4 to 8 Range 

Linear Theory - Supersonic linear theory is developed from 
small perturbation considerations. Several of the key assump- 
tions leading to the linearized potential flow equation effec- 
tively limit the Mach range and vehicle geometry for which the 
theory applies. These assumptions are: (1) that the thickness 

and camber slopes are small, (2) that the products of the pertur- 
bation velocities are negligible compared to linear terms in the 
perturbation quantities, and (3) that the flow is irrotational 
or isentropic. Assumption (1) places limitations on the vehicle 
geometry. Supersonically, it means that no slope on the vehicle 
can be larger than the Mach angle. This restriction makes it 
very difficult to model real configurations at hypersonic Mach 
numbers. The second assumption is necessary to linearize the 
equations, but it also places an upper limit on the Mach range, 
since the neglected terms are all multiplied by M^ . The third 
assumption, that the flow is isentropic, implies that shock 
waves can be approximated by Mach waves. Since the Mach angle 
is, in general, less than the shock angle, the effect is two- 
fold. The strength of the wave is incorrectly computed, and the 
region of influence, i.e., within the Mach cone, is underpredic- 
ted. The relative error in the region of influence increases 
rapidly with Mach number, as shown for a 10° cone in Figure 1. 
Since the effectiveness of wing camber in reducing drag-due-to- 
lift depends on the region influenced by each element, this is 
one of the reasons why unmodified linear theory applied in the 
inverse design mode shows little or no payoff above Mach 5. 


The supersonic linear theory formulations of References 1 
to 3, as applied to wings, rely on solutions which uncouple the 
camber and thickness terms. These solutions neglect the contri- 
butions of the thickness slopes to the local lifting pressures. 
The present study shows that these terms are not negligible 
above Mach 4, and that they must be included, along with non- 
linear camber terms, for appropriate aero predictions in the 
Mach 4 to 8 range. 


2 




FIGURE 1 

COMPARISON OF REGION OF INFLUENCE AS DEFINED BY MACH ANGLE 
(LINEAR THEORY) AND SHOCK ANGLE ON 10° CONE 


The result is that linear theory methods strictly apply 
below about Mach 3 for configurations that may be characterized 
as sharp, thin, with small camber, and at low angle of attack. 

In addition, various practical problems arise in the trans- 
lation of the theory into operational computer programs, even in 
the applicable Mach range. These problems relate to the trans- 
lation of the geometry of a complex, 3D vehicle into an appro- 
priate linear theory analog for the computer program. 

Impact Theory - Newtonian impact theory is based on the 
assumption that the Mach number is so high (M «=) that the shock 
is nearly coincident with the surface, and that the pressure on 
the surface is, therefore, only a function of the local surface 
slope. The implicit assumption is that the Mach angle, and 
therefore the region of influence from a point, is zero. The 
result is that impact theory is really two-dimensional. Given a 
surface element, no account is taken of the effect of adjacent 
elements, or of the element geometry (sweep, aspect ratio, taper 
ratio, or tip effects). Thus, the three dimensionality of the 
vehicle planform and the interference effects are not properly 
assessed unless the free stream Mach number is very high. 
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Attempts have been made to extend impact theory to lower 
Mach numbers by introducing a Mach number -dependence, as 
well as surface slope, into the pressure coefficient relation. 

The result is that there are now some 27 different pressure 
options to choose from in the Gentry program. While the careful 
choice of the pressure options applied to complex configurations 
in the Mach 4 to 8 range generally gives good force and moment 
predictions, as a rule, the local pressure distributions are not 
well predicted. For this reason, attempts to optimize high 
speed configurations using impact theory methods have not been 
very successful. 

Combined Impact Theory/Linear Theory Approach 

Impact theory is conceptually and computationally simple, 
and permits a geometric representation of considerable detail. 

The method correctly accounts for the nonlinear aerodynamics 
which characterize the higher Mach range. However, the approach 
lacks the mechanism for properly assessing both the effect of 
the local element geometry and the influence of one surface ele- 
ment on another. This mechanism exists in the aerodynamic 
influence coefficient (AIC) matrix of a number of finite element, 
linear theory analyses. 

The aero influence coefficients, obtained from the solutions 
to the linearized differential equations governing the supersonic 
flow field, give the effect of a disturbance on the i-th element 
and on all other elements of the configuration. Using a chord 
plane distribution of vortex elements, as in the analysis of 
Reference 1, the basic linear theory equation can be written 

[a ij ] {AC pj } - {a ± } (1) 


In the above matrix equation, a. . is the coefficient for the 
influence of panel i on panel j^ACp- is the pressure difference 
(vortex strength) across the j-th panel, and is the linearized 
boundary condition on panel i. In this case, af is the sum of 
the angle of attack, the local camber slope, and the interference 
onset flow from all other sources, i.e.. 


a . = a — a + n . 

1 C . 1 

1 


( 2 ) 


Solving equation (1) for the lifting pressures, which requires 
inverting the AIC matrix. 


or 




AC = I aT^ a. 
Pj i H a 


(3) 
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( 4 ) 


Now, the lifting pressure on a 2-D flat plate is given 



where 6 is /m 2 -1 and the asterisk denotes the 2-D non-interference 
pressure. Substituting aj_ from equation (4) into (3), we see 
that the role of the influence coefficients is to correct the 
2-D non-interference pressures for three dimensionality and in- 
terference effects, i.e., 

8 —1 * 

AC = ~ Y. a. 7 AC (5) 

Pj 4 . 13 p. 

Equation (5) is the basic equation for the combination of 
linear theory and impact theory. Using appropriate impact pres- 
sure options, the Gentry program is used to calculate the 2-D, 
noninterference pressures over the configuration. The aero- 
dynamic influence coefficients from the Woodward program are 
then used to modify the upper and lower surface impact pressures 
using equation (5) . The modified pressures are then integrated 
over the Gentry geometry model for the forces and moments. 

The appropriate impact pressure options to be used in this 
analysis are those which give the correct local compression and 
expansion pressures on an inclined 2-D flat plate. These are 
the oblique shock (or tangent wedge) relation for compression 
pressures, and the Prandtl-Meyer relation for expansion pres- 
sures. For blunted configurations, a pressure option which ex- 
hibits Newtonian flow characteristics at deflection angles near 
90°, and which fairs smoothly into the oblique shock solution 
at low angles is needed. The following pressure relation exhibits 
the desired characteristics. 


r = 9 sin9 sinS 
p cos ( e — 5 ) 

where (6) 

6 = p + <5 - p sin 6 


and 6 is the local flow deflection angle and u is the freestream 
Mach angle. 

VALIDATION OF COMBINED IMPACT THEORY/LINEAR THEORY METHOD 
70° Swept Flat Plate at Mach 6.0 

The combined Gentry/Woodward analysis was used to calculate 
the pressures on a 70° swept flat plate at Mach 6 and 8° angle 
of attack for comparison with unpublished NASA data. The results 
of the straightforward application of the approach are compared 
with the data in Figure 2. Also shown in Figure 2 are the 
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Woodward -a lone (pure linear theory) results, and the Gentry-alone 
(pure impact theory) results. The comparisons show that the 
spanwise pressure distribution calculated using the combined 
Gentry/Woodward approach is significantly better than either the 
basic linear or impact theories. 
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FIGURE 2 

PRESSURE COMPARISONS ON 70° SWEPT FLAT PLATE 

Mach 6.0 a - 8° 


Although the straightforward application of the combined 
approach gives the proper trend and approximately the correct 
pressure level on the flat plate, the spanwise pressure distribu- 
tion is shifted too far inboard. This is because the influence 
coefficients calculated in the Woodward program are based, in 
part, on the assumption that the region of influence for a given 
element is contained within a Mach cone corresponding to the 
freestream Mach number of 6.0. Actually, the average local Mach 
number, corresponding to the average local pressure on the flat 
plate, is close to 5.0, and the region of influence based on the 
average local Mach number is about 20 percent greater than given 
by the freestream Mach cone. The influence coefficients for 
this case were recomputed where a spanwise scale factor of 1.2 
was used to decrease the span for the linear theory calculation. 
This spanwise scale factor effectively increases the region of 
influence to correspond to the average local Mach number on the 
planform. Figure 3 shows that the use of a spanwise scale fac- 
tor to correct the region of influence gives improved results. 
These results support the previous assertion that unmodified 
linear theory applied at high Mach numbers will tend to under- 
predict the benefits due to wing camber because the region of 
influence is underpredicted. 


6 



0.10 



FIGURE 3 

EFFECT OF SCALE FACTOR ON GENTRY/WOODWARD SOLUTION 
FOR 70° SWEPT FLAT PLATE 

Mach 6.0 a= 8° 


76° Swept Wing Comparisons 

Experimental pressure data is given in Reference 6 on a 76° 
swept wing with thickness and camber ( c L^ es ^g n = 0.0) for several 

Mach numbers up to 4.6. The combined Gentry /Woodward program was 
used to calculate the pressures on this wing at Mach 2.3, 3.5, 
and 4.6. The results at Mach 4.6 are compared with the data for 
three angles of attack at six different spanwise stations in 
Figures 4(a) through 4(f). Also shown are the pressure distribu- 
tions calculated using both Gentry-alone and Woodward-alone. The 
Woodward-alone results include the use of a vacuum pressure 
limit calculation. 

At the centerline. Figure 4(a), the pressures calculated 
with the combined Gentry /Woodward analysis give excellent compari- 
sons with the data, on both the upper and lower surfaces. As 
might be expected, the comparisons are best for the low angle of 
attack (5.5°), but even at the highest angle of attack (20.5°), 
which is nearly twice the freestream Mach angle, the results are 
very good. As one moves outboard, Figures 4(b) to 4(f), the 
combined theory begins to over-predict the pressures at the 
leading edge, and this over-prediction grows progressively worse 
toward the tip. At the same time, the Gentry-alone pressures 
compare better and better with the data as one moves outboard. 
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FIGURE 4 (a) 

PRESSURE DISTRIBUTIONS ON 76° SWEPT WING 

Mach 4.6 Span Station (y/b/2) = 0.0 


a = 5.56° a =10.57° a = 20.57° 
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FIGURE 4 (b) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 4.6 Span Station (y/b/2) = 0.2 
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FIGURE 4 (e) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 4.6 Span Station (y/b/2) = 0.8 


a = 5.56° a = 10.57° a = 20.57° 



FIGURE 4 (f) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 4.6 Span Station (y/b/2) = 0.95 
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The over-predictions in the region of the leading edge re- 
sult from the linear theory influence coefficients. Comparisons 
with the closed form linear theory solutions for triangular flat 
plates, given in Reference 7, show that the influence coefficients 
from the Woodward program are correctly calculated. The indi- 
cation is that the linear theory solution is inappropriate in 
the region of the leading edge when the planform has thickness 
and the leading edge is swept near or aft of the Mach cone. It 
is believed that this problem arises from the linearization of 
the basic supersonic flow equation. This problem can be expec- 
ted to persist until an improved method of calculating the in- 
fluence coefficients becomes available. 

Inspection of the Gentry /Woodward and Gentry-alone compari- 
sons in Figure 4 gives rise to the following empirical rule for 
correcting the leading edge over-predictions: In the region 

where the Gentry /Woodward compression pressure is greater than 
the Gentry-alone pressure, use the Gentry solution, otherwise 
use the Gentry/Woodward solution. Comparison with the data 
shows that application of this rule gives reasonably good re- 
sults at all wing stations. Figure 5 shows that, even though 
the local pressures are not uniformly well predicted at all 
points, the method gives excellent predictions of the pressure 
differences at all wing stations. In addition, this simple rule 
has the advantages of being locally determined and easily pro- 
grammed . 
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FIGURE 5 (a) 

PRESSURE DIFFERENCES ON 76° SWEPT WING USING CURRENT METHOD 

Mach 4.6 a =5.56° 
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FIGURE 5 (b) 

PRESSURE DIFFERENCES ON 76° SWEPT WING USING CURRENT METHOD 

Mach 4.6 a = 5.56° 


Application of the above rule also provides a method for 
correcting the uninverted aero influence coefficient matrix ele- 
ments, which are used in the inverse design method. Using 
equation (5), the Gentry/Woodward pressure on the j-th element 
is computed from 



* — i 

where Cp^ is the impact pressure on the l-th element, and a^j is 
the influence coefficient which gives the contribution of the 
i-th element on the j-th panel. If the Gentry/Woodward pressure 
on the j-th element is greater than the Gentry-alone pressure, 
i.e., Cp.>Cp., we assume that we can scale the influence coeffi- 
cients aJfect^ing the j-th panel so that 



E 

i 



— C 

SCALED 4 p. 


( 8 ) 


Assuming that the same scale factor is applied to each coeffi- 
cient affecting the j-th panel, this means that the j-th row 
in the inverted matrix is multiplied by a constant, and there- 
fore, the j-th column in the original uninverted matrix is 
divided by that constant. Note that only the influence coeffi- 
cients affecting the j-th panel are modified. 
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Additional Comparisons on the 76° Swept Wing 

The combined Gentry/Woodward analysis was also used to cal- 
culate the pressure distributions on the 76° swept wing at Mach 
2.3 and 3.5. The comparisons with the data at Mach 3.5 are given 
in Figure 6, and the Mach 2.3 results are given in Figure 7. In 
general, the comparisons are similar to the Mach 4.6 results, but 
the quality of the local pressure predictions decreases with de- 
creasing Mach number. Since the combined solution at the lower Mach 
number (2.3) is very close to the pure linear theory solution, as 
shown in Figure 8, the poorer pressure comparisons are attributed 
to the linear theory formulation for subsonic leading edges. Never- 
theless, the use of the linear theory influence coefficients 
to correct the impact theory pressures allows the Gentry program 
to be used to calculate the forces and moments at these low 
supersonic Mach numbers, as illustrated in Figure 9. 


a = 5 . 73 ° 


a « 9 . 74 ° 


a= 19 . 74 ° 



FIGURE 6 (a) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 3.5 Span Station (y/b/2) = 0.0 
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FIGURE 6 (d) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 3.5 Span Station (y/b/2) = 0.6 
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FIGURE 6 (e) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 3.5 Span Station (y/b/2) = 0.8 




a = 5.73 a = 9.74 a = 19.94 



FIGURE 6 (f) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 3.5 Span Station (y/b/2) = 0.95 



FIGURE 7 (a) 

PRESSURE DISTRIBUTIONS ON 76° SWEPT WING 

Mach 2.3 Span Station (y/b/2) = 0 


16 



a - 5 . 94 ° a - 9 . 94 ° a - 19 . 94 ° 



FIGURE 7 (b) 

PRESSURE DISTRIBUTIONS ON 76° SWEPT WING 

Mach 2.3 Span Station (y/b/2) = 0.2 


c. - 5.94° a = 9.94° a = 19 94° 



FIGURE 7 (c) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 2.3 Span Station (y/b/2) = 0.4 
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FIGURE 7 (d) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 2.3 Span Station (y/b/2) = 0.6 
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FIGURE 7 (e) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 2.3 Span Station (y/b/2) = 0.8 




FIGURE 7 (f) 

PRESSURE DISTRIBUTION ON 76° SWEPT WING 

Mach 2.3 Span Station (y/b/2) = 0.95 



y/b/2 = 0.95 
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FIGURE 8 

COMPARISON OF LINEAR THEORY, IMPACT THEORY AND COMBINED APPROACH 

FOR PRESSURES ON 76° SWEPT WING 

Mach 2.3 a = 5.94 
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FIGURE 9 

LIFT AND MOMENT ON 76° SWEPT WING AT MACH 2.3 


Ames All-Body Models 

The combined impact theory/linear analysis was applied to 
the NASA Ames All-Body models illustrated in Figure 10. The 
Gentry /Woodward program was used to calculate the lift and moment 
on the elliptical cone model (B 5 ) at Mach numbers 2.0 and 5.37. 
The predictions are compared with the data from Reference 8 in 
Figures 11 and 12 and show that the combined analysis can be 
applied across the Mach range to obtain results equal to, or 
better than, either linear theory or impact theory used alone. 

The lift comparisons at both Mach numbers indicate that there is 
a thickness contribution to the lift which is not predicted by 
linear theory. The importance of these thickness contributions, 
particularly as they relate to an inverse design method in the 
Mach 4 to 8 range, are discussed below. 
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BODY B-| - ELLIPTICAL FOREBODY/CLOSED AFTERBODY 


B = 0.4826 m 

s = 0.06241 m 2 
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BODY B 5 - ELLIPTICAL CONE (SAME AS 



B-| FOREBODY) 
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a = 0.0216 m 

b = 0.0862 m 


XR£p = 0.177 m 
c = 0.215 m 


NASA AMES ALL-BODY MODELS B-] AND B 5 



FIGURE 11 

LIFT AND MOMENT PREDICTIONS ON AMES ALL-BODY MODEL B 

Mach 2.0 




FIGURE 12 

LIFT AND MOMENT PREDICTIONS ON AMES ALL-BODY MODEL B 5 

Mach 5.37 


Closed Form Solutions for Lift, Drag, and Moment 

Starting with the basic relation for combining linear theory 
and impact theory, the impact pressure differences obtained from 
the small angle expansions of the oblique shock and Prandtl-Meyer 
relations wete substituted into Equation (5) . The equation was 
then expanded further to give the lifting pressure in terms of 
the influence coefficients, the slope of the camber line measured 
from the freestream, and the local thickness slope. The result- 
ing expression for the lifting pressure was then summed over the 
planform for the lift, moment, and drag-due-to-lif t . With the 
exception of higher order terms in local camber slope, the re- 
sulting summations could be identified as the linear theory con- 
tributions to the lift, drag, and moment, or in the case of the 
thickness terms, using fairly general approximations, could be 
replaced by integrable expressions. From these expansions, then, 
closed form expressions were derived for the lift, drag, and 
moment for uncambered bodies. These closed form expressions, 
given in Figure 13, require only the lift and moment curve slopes 
as calculated from linear theory, the angle of attack, and the 
base area and volume of the configuration. 
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FIGURE 13 


CLOSED FORM EXPRESSIONS FOR LIFT, MOMENT AND DRAG-DUE-TO-LIFT 
FOR UNCAMBERED CONFIGURATIONS 


The closed form expressions were used to calculate the char- 
acteristics at Mach 5.37 and 7.38 on the Ames All-Body models 
illustrated in Figure 10. The results are compared with the 
data in Figures 14 and 15. These results show that if the con- 
figuration has a base area, then there is a thickness contribu- 
tion to lift curve slope which is proportional to the base area. 
Further, even if the configuration is closed at the base (zero 
base area) , there is a thickness contribution to the moment curve 
slope (Cm ) which is proportional to the volume. These results 
demonstrate that, in general, the thickness terms must be in- 
cluded in the inverse design procedure to ensure that the opti- 
mum camber distributions correspond to the actual lift and moment 
constraints . 
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FIGURE 14 
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FIGURE 15 

CLOSED FORM PREDICTIONS - AMES BODIES B-| AND B 5 
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INVERSE DESIGN PROCEDURE 


The use of linear theory influence coefficients to correct 
the impact pressure distributions provides the basis for the 
use of the method in the inverse design mode. Three things must 
be taken into account for application in the Mach 4 to 8 range: 

1) the influence coefficients in the region of the leading edge 
must be corrected to prevent the large leading edge overpredictions 
in the lifting pressures. 2) The thickness contributions to the 
lifting pressures must be included in the analysis. And, 3) the 
region of influence must be corrected for large values of 8 cot A, 
where A is the leading edge sweep. 

The method for scaling the influence coefficients in the 
region of the leading edge by comparing the linear theory pres- 
sure to the impact pressures was previously discussed. Also, 
the use of a spanwise scale factor to correct the linear theory 
region of influence was demonstrated for the 70° swept flat 
plate at Mach 6.0. The remaining task is to show that the thick- 
ness contributions to the lifting pressures at high Mach numbers 
can be introduced into the linear theory inverse design procedure. 

The small angle expansions of the oblique shock and 
Prandtl-Meyer relations can be used to give an expression for 
the impact theory lifting pressure in terms of camber line slope, 
thickness slope, and Mach number. This expression is 

AC* = £ a. (1 + 1 . 26 e • + .63 3/2 e. 2 ) + .83 1/2 a. 3 (9) 

p . p i l l l 


where a . is the slope of the camber line measured from the free- 
stream, ej_ is the local thickness slope as calculated from the 
Gentry program for a ^ = 0, and 6 is /m 2 - 1 . Substituting equa- 
tion (9) into the basic equation, equation (5), and neglecting 
the cubic term in a-^, the lifting pressure is 



a . . a ■ b • 
i j ii 


( 10 ) 


where bi is the thickness term in equation (9) . With this ex- 
pression, the optimization equations derived by Woodward in 
Reference (1) for minimum drag can be rederived, to include the 
thickness terms. The resulting set of linear equations, for 
the optimum loadings is given, in matrix form, in Figure 16. 
Solving this set of equations for the ( ACpj ) OPTIMUM* the corres- 
ponding optimum camber distributions are calculated from 


a 


i 


1 _ 

b. 


Z. 

j 


a . . 


13 



OPT 


(ID 


Note that for zero thickness, b^ -»■ 1, and the optimization equa- 
tions reduce to the same set derived by Woodward. 
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OPTIMIZATION MATRIX - Solution gives loadings for minimum drag 
corresponding to lift and moment constraints 


s l a ll s l a ll 
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OPTIMUM CAMBER - calculated from optimum loadings 


where 


a . . 
iD 

b. 

1 


A C p . 

1 

x-x . 

1 

A l' A 2 
L, M 


a i = b7 E a ij ( AC p . } OPT 
ID D 

= area of element i 

= influence coefficient (element i on element j) 

3/2 2 

= 1 + 1.2(3e^ + e . (thickness factor from 

Gentry Program) 

= pressure difference on element i 

= distance from moment ref. to centroid of element i 
= langrangian multipliers 
= design lift and moment constraints 


L=C t S M=C m Sc 

L M 

FIGURE 16 

OPTIMIZATION MATRIX FOR DESIGN PROCEDURE 
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For large values of camber slope, and/or for very high Mach 
numbers, the cubic term in camber slope may not be negligible. 
For this case, an iterative procedure may have to be applied to 
obtain the actual camber slope. If iteration is necessary, the 
initial values of the bj^ factor would be adjusted using the pre- 
viously calculated camber slope, so that 


(b i ) iteration ^i^ + ' 20 


3/2 


) 


previous 


where (bi) 0 is the thickness term neglecting the higher order 
camber term. 

The procedure for applying the combined impact theory/linear 
theory approach is outlined in Figure 17. 
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STEPS IN DESIGN PROCEDURE 
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CONCLUDING REMARKS 


The results of this study demonstrate that the combination 
of linear theory and impact theory is a feasible and attractive 
approach for improved calculations of the aerodynamic character- 
istics of configurations at high speeds. The application of the 
combined approach to several wing alone cases shows that the 
method not only gives improved predictions for local pressures 
and loadings, but that it also allows the basic Gentry program to 
be used at low supersonic Mach numbers to predict forces and 
moments. Additional benefits of the combined approach are 1) it 
reduces the number of impact pressure coefficient options required 
to just three basic options, 2) it provides a method for correc- 
ting the linear theory influence coefficients to remove the over- 
predictions in lifting pressures near the leading edge when the 
planform is swept near the Mach cone, and 3) it provides the basis 
for an inverse design procedure applicable in the Mach 4 to 8 
range. 

In the course of the study, the basic relations for the com- 
bination of linear theory and impact theory were used to derive 
closed form expressions for the lift, moment, and drag-due-to- 
lift of uncambered bodies in terms of the linear theory slopes, 
the angle of attack, and the base area and volume of the configur- 
ation. These expressions should be useful for the preliminary 
estimation of aerodynamic characteristics, including thickness 
and non-linear lifting effects, across the Mach range. 

Although the procedure has not been applied to wing-body con- 
figurations, the results obtained in the present study on several 
highly swept wings (up to 76°) indicate that extension to wing- 
bodies is feasible. Two linear theory methods are available for 
calculating the influence coefficients for the body and body flow 
field contributions. These are the use of chord plane vortex 
panels to represent the body, in the manner similar to that 
described in reference 2, or the use of surface singularities as 
was done in reference 3 . 
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APPENDIX 


CLOSED FORM SOLUTIONS FOR LIFT, DRAG, AND MOMENT 
FOR UNCAMBERED CONFIGURATIONS 


Basic Equation - The basic equation for the combination of 
linear theory and impact theory is 



(Al) 


where ACpj is the lifting pressure on element j, a^j - is the in- 
fluence coefficient for panel i on element j, ACpi is th e 2-D , 
noninterference lifting pressure on element i, and 6 = /m^-1 . 


Impact Pressures - The 2-D, noninterference pressures are ob- 
tained from small angle expansions of the oblique shock and 
Prandtl-Meyer relations. 


C* = 1 6. + 1.2 6? + .4 ✓& 6? (A2) 

p^ e i i i 


in the above expression 6^ is the local surface slope, and 


( 6 i ) LOWER 
* 6 i^ UPPER 


e . + a . 
i l 


e . - a . 


l i 


(A3) 


where is the upper surface thickness slope measured from the 
camber line, and aj_ is the slope of the camber line measured 
from the freestream. If a c ^ is the slope of the camber line 
(dz c /dx) measured from the reference chord plane, then 


(A4) 


With the above expressions, the 2-D, noninterference lifting 
pressures, and the thickness pressures are, respectively 


* * 
AC = (C ) 


* 

,-(C ) 


p . ' LOWER p ■ UPPER 


= a . + 4.8a.e. + . 8/3 (3e? + a?)a. 


1 1 


and 


(A5) 


V 'lower + < C p . ’upper = f E 1 + 2 - 4 <“i + e i> + - 8/B U i + 
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Substituting the impact pressure relations (A5) into the basic 
equation (Al) , the combined approaches gives 


AC = E a. 1 a . + .23 3 ^ 2 E a . \ a? + 1.23 E aT ^ e . ct . 

Pj j, ID i i ID i ± ID i i 

. a 3/2 _ -1 2 

+.63 E a . . e . a. 

. r D r r 

(C p i ) LOWER + (C p i ) UPPER ? a ij e i + ,6B f a ij e i 


3/2 -13 -12 3/2 -1 ? 

+ .26 7 E a., e 7 + .63 E a . 7 a + .66 7 E a . . e.a. 


ID i 


ID i 


13 ii 


Lift, Drag, and Moment - The lift, drag, and moment coefficients 
are obtained by summing the pressures over the planform areas. 
Thus 


C L S = ' S j AC p. - E S j [(C p. } UPPER + (C P . } LOWER 1 
D D D D ^D 


'p. “j + E S D [(C P .> UPPER 
*3 D D 


+ ^p.^ LOWER-* E j (A7) 


C M S C = E S . AC (x - x . ) 

M j p . 3 

D D 

where S is the planform area, Sj is the area of panel j , C is 
the reference chord, and (x - x j ) is the distance from the 
moment reference center to the centroid of panel j . 

Substituting the pressure relations from (A6) into the force 
and moment relations (A7) we have 
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C T S = Z S. E a.-!- a. + .2g 3//2 E S. E a . "I - a? 

L j ^ i ^ 1 j ^ i ^ 1 

+ 1 . 2 g E S. E aT^ a.e. + .6g 3 / 2 E S. E aT* e? a. 

j 3 i XD i i j '3 i 13 i x 


-1 


HIGHER 


-ES. e. a. E a.t e. + (ORDER ) 


3 3 3 


13 i 


TERMS 


1 3/2 -1 3 

C_. S = E S. a. E a., a. + .2g ' E S. a. E a., a. 

D j 3 3 ± 13 i j 3 3 i 13 i 

3/2 -1 2 -1 

+ . 6 g ' E S. e. E a., e.a. + 1.2g E S. a. E a . 7 a.e. 

j 3 3 L 13 i i j 3 U i 13 11 

-1 2 3/2 -1 2 

+ .6g E S. e. E a.j a. + • 6 g E S. a. E a., e. a. (Ag) 

+ E S. e. E aT'!' e. + .6g E S. e. E aT'!' e 2 

j 3 3 i 13 1 j 3 3 i 13 1 


+ .20 


3//2 E S . e . E aT * e? 
j 3 3 i 13 1 


C.. S C = E S. (x - x.) E a. 1 a. + .2g 3 / 2 E S . (x 
M j 3 3 L 13 1 j 3 


v _ -1 3 

x . ) E a . . a . 

3 i 13 1 


+ 1.20 E S. (x - x.) E a.* a.e. 

3 3 ± 13 1 1 


3 

3/2 


+ .60 ' E S . (x - x.) E a., a. e. 

■ 3 3 • 13 1 1 

j J j -l -1 


Identification of Terms - The first summation in each of the 
above equations is the linear theory contribution for a cambered 
flat plate. Thus 
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If the plate is uncambered, then = a and these terms become 
(LIFT) 


E S. I a. . a. -*■ a E S . E a . . — a (C T ) T ~ S 

j : i 1 ] 1 j 3 i 13 L a LT 


(A10) 

UNCAMBERED 


(DRAG) ES.a.Ea.^a.-*ct 2 ES.E -1 2 , e 

it x i . a. . = a C ) S 

i JJ l J J J i li L LT 

J J J a 

(MOMENT) E S.(x-x.) E aT 1 a . -> a E S • (x-x . ) E aT 1 = ct (C„ ) Tm S C 
j 3 3 i II i j J J ± 13 M a LT 


The second term in each of the expressions in (A8) represent 
the non-linear contributions of angle of attack and camber slope. 
For an uncambered configuration, where a. = a, these terms be- 
come 
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(C M ^ LT S 

a 
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In 

the drag 

expression 

(A8), the 

last three 

terms give the 

drag due 

to thickness (or 

wave drag) . 

We will simply replace 


these terms by Cd-thiCK Si 


Uncambered Expressions (a-j = a) - Using the above expressions, 
the equations (A8) are rewritten 
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C T S = (C T ) Trp S a (1 + . 20 3/2 a 2 ) + 1.26 a E S. E a . ^ e. 
L L a LT j ^ i 13 i 

3/2 -1 2 -1 

+ .66 aES. Ea.. e. — a E S . e. Ea.. e. 

j 3 i 13 i j 3 3 ± 13 i 


(C D- C D, 


THICK 


< C L >LT S “ 2 

(1 + .26 3/2 a 2 ) 
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+ .66 3//2 a 2 [ E S. E a~^ 

j 3 i 13 1 

j 3 i 13 


-1 2 -1 
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3 3 . 13 l 3 3 . 13 

3 J J 1 J 3 J J 1 J 

C M S C = (C M ) LT S C a(l + . 2 6 3/2 a 2 ) 

a 

+ 1.26 a E S. (x-x.) E a7 3 e. 

j 3 3 ± 13 1 

+ .66 3/2 aES. (x-x.) E aT^ e 2 
j 3 3 j 13 1 

Thickness Terms - The summations in equation (A12) which are 
linear m thickness slope can be evaluated under certain condi- 
tions. For instance, if the thickness slope is constant over 
the configuration: 


ES. E a., e. = e ES. E a.^ = e S (C T ) T m 
■ 3-131 0-3-13 o L L.T. 

j 1 J 3 J 1 J a 


-BASE 


(C L ^ LT 
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At very high Mach numbers (M °°) , there is no influence 
of one element on another. For this case, the influence coeffi- 
cients are 

a . } = 0 for i ^ j, and a . ^ = % for i = j 
13 J 136 
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where the linear theory lift curve slope approaches 4/3 at hyper- 
sonic Mach numbers. 

The same result is obtained for an unswept wing at all su- 
personic Mach numbers. 

From the above considerations, we will approximate the in- 
fluence coefficients in the summations involving thickness by 

aT^ = 0 for i / j , and aT 3 : = (C ) T for i = j 
3 3 a 

Then the double summations in equation (A12) can be approximated 
by single summations. 


C L S = (C ) S a(l + .2B 3/2 a 2 ) + 1.23 a (C L ) L>T> E S. e. 

a a * * j J J 

+ ( . 63 3/2 - 1) a(C T ) T m E S. e • 2 

L a LT j 3 3 

(C -C ) s = (C_ )__ s a 2 ( 1 + .23 3/2 a 2 ) 

D D THICK L a LT 


(A13 ) 


+ 1.26 « 2 (C ) I S c + .66 a 2 (C ) I S e 

a j a j J J 

3/2 2 2 

+ ^ a < C L } LT E S j £ j 
a j j j 

C S C = (C M )__ S C a(l + .23 3/2 a 2 ) + 1.23 a (C_ ) TT E S. (x-x ) e. 
•M M a LT L a LI . J J J 

+ .63 3/2 a (C )__ E S (x-x.) e. 2 

L a LT j ] 3 3 

Taking the thickness slope to be (dz/dx) and replacing Sj by 
dxdy, we write 


b/2 TE 


dz , 


b/2 


BASE 


E S. e . - / / <g£) dxdy = / Z dy - 

j 3 3 -b/2 0 d -b/2 


Similarly 


rr A BASE 


E S . (x-x . ) e . = X E S . e.-ES. X. e . — X 0 
j 3 3 J j J J j J J 3 2 


- E S . X . e . 

j 3 3 3 


34 



The summation £ Sj Xj is replaced by an integral and integrated 

by parts, where ej is replaced by f' and f is the equation for 
the upper surface. Thus 


HDD 


»/2 v 

—b/2 XtE EtE dY " 2 

where f T „ is assumed to be zero and V is the volume. 

Li-Cj 

The remaining integral is X__ where X mTP is the dis- 

I Ji dAoiIj llli 
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Z Sj (X-Xj) Ej = ( X “X TE ) 


BASE , V 
2 2 


We have assumed that e j is the thickness slope in the x-z 
plane. Actually, ej should be the flow tangency angle. For 
instance, for a cone of semi-vertex angle e c , 

£ S . E . = 8 £ S . = 0 S = 2 (^) = F(^SE) 

in c.i c it 2 2 

J J 

Similar considerations for other simple shapes indicate the 
deviation from our basic 2-D assumption can be corrected by in- 
troducing a factor F which is 

F = s /<! E wet> 

where S is the planform area, and S wet is the wetted area of 
the configuration. 

With the above considerations, the lift, moment, and drag 
equations are 


C L = <C L >LT “ U + ‘ 2eV 2 “ 2 » + •« “ (C L >LT F 
a a 

+ (.6g 3/2 - 1) a (C L ) LT (e 2 ) 

“a 


(A14 ) 
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(Cd " Cd THICK ) = C L a + - 3e “ 2 (C L } LT F 

a 

+ ( . 60 3/2 + 1) a 2 (C L ) LT (e 2 ) 

a __ 

. ..3/2 2. , ,. ,. , _ , X TE A EASE , V . 

C M - < C M LT 0,(1 + - 26 a 1 + a C L LT F <— 5“ + ST> 

a a C 


+ .6g 3/2 a(C L ) LT 


Z s . (X-X . ) e . ' 
2 2 2 


and 



2 

Comparisons of the e terms against the other terms in- 
dicate that for most configurations these terms can be neglected. 
However, given the details of a configuration, the contribu- 
tions of these terms can be estimated. 
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